arXiv:1503.03777vl [cond-mat.mtrl-sci] 12 Mar 2015 


Control of Dzyaloshinskii-Moriya interaction in Mni_j;Fea;Ge: a first-principles study 


Takashi Koretsune/ Naoto Nagaosa/’^ and Ryotaro Arita^ 

^RIKEN Center for Emergent Matter Science (OEMS), Wako, Saitama 351-0198, Japan 
^Department of Applied Physics, University of Tokyo, Hongo, Tokyo 113-8656, Japan 

(Dated: March 13, 2015) 

Motivated by the recent experiment on the size and helicity control of skyrmions in Mni_ 3 ;Fea;Ge [K. Shi- 
bata et ah, Nature Nanotechnology 8, 732 (2013)], we study how the Dzyaloshinskii-Moriya (DM) interaction 
changes its size and sign in metallic helimagnets. By means of first-principles calculations, we successfully 
reproduce the non-trivial sign change of the DM interaction observed in the experiment. While the DM inter¬ 
action sensitively depends on the carrier density or the detail of the electronic structure such as the size of the 
exchange splitting, its behavior can be systematically understood in terms of the distribution of anticrossing 
points in the band structure. By following this guiding principle, we can even induce gigantic anisotropy in 
the DM interaction by applying a strain to the system. These results pave the new way for skyrmion crystal 
engineering in metallic helimagnets. 


A skyrmion is a topologically protected nano-size spin tex¬ 
ture found in several magnets^— . Due to its unusual spin struc¬ 
ture, many intriguing behaviors such as topological Hall ef¬ 
fects, current-driven motion, and multiferroic behavior have 
been observed^ Although there is a huge potential to design 
novel functional materials by exploiting these unique elec¬ 
tromagnetic properties, skyrmion engineering or skyrmion- 
crystal engineering is yet to be established. Here, the key is¬ 
sue is how to manipulate the size and helicity of skyrmions. 
Regarding this problem, a recent experiment for the represen¬ 
tative skyrmion system Mni_2:Fea;Ge has shown that we can 
tune the skyrmion size and helicity by changing the carrier 
density^. 

The Hamiltonian which determines the nature of skyrmions 
is 


H = 



{VMf -b DM • (V X M) , 


( 1 ) 


where M is the magnetization per volume, D is the 
Dzyaloshinskii-Moriya (DM) interaction coefficient and J is 
the ferromagnetic exchange coupling, respectively. For mate¬ 
rials design of skyrmion crystals, we need to know the precise 
value of D/ J. However, non-empirical evaluation of these pa¬ 
rameters in the classical continuum model has been a difficult 
challenge, since it requires an elaborate multi-scale approach 
spanning the quantum to classical regime. For the insulating 
skyrmion system Cu20Se03^, there is a work in which the 
spin Hamiltonian ( 1 ) was derived from first principlesiS. By 
comparing the total energy of various magnetic states, they de¬ 
termined the value of D/J, and succeeded in reproducing the 
experimentally measured skyrmion size. However, the guid¬ 
ing principle to control the values of D/J is yet to be ob¬ 
tained. 

On the other hand, for metallic systems, the situation is dif¬ 
ferent. There are several studies for the estimate of D^—. 
Among them, recently, one of the present authors (NN) and 
his collaborators have shown that the DM interaction in the 
two-band model drastically changes when the band anticross¬ 
ing point resides near the Fermi level^li. The story is anal¬ 
ogous to that of anomalous Hall conductivity, in which the 
band anticrossings act as magnetic monopoles in momentum 


spaced. In Ref. [l^ a Berry phase expression for the DM in¬ 
teraction has also been formulated. These studies stimulate 
us to explore a fascinating possibility of controlling the DM 
interaction in metallic systems by manipulating the electronic 
structure. Indeed, the fact that not only the size but also the 
helicity of skyrmions in Mni_a;Fe^Ge changes as a function 
of X indicates that we have a good chance to control the value 
ofD. 

In this paper, we show a quantitative analysis of the DM in¬ 
teraction in the metallic helimagnet, Mni_a;Fe^Ge, based on 
ab initio density-functional theory (DFT) calculation. From 
the obtained band structure, we evaluate the off-diagonal 
spin susceptibility which is a direct measure of the DM 
interaction. We find that the sign change of D observed in 
the experiment for Mni_a;Fea;Ge is successfully reproduced. 
The carrier-density dependence of D can be systematically 
understood in terms of the distribution of band anti-crossing 
points in the electronic structure. We demonstrate that the 
sign and the size of D can be controlled as a function of the 
carrier density or the size of the exchange splitting. There is 
also an interesting possibility to induce gigantic anisotropy in 
D by applying a strain to the system. 

Results 

DM interaction in the continuum model. Let us first look 
at the second term in the Hamiltonian ([T]i. This indicates that 
g-linear term in the spin susceptibility, should be pro¬ 
portional to the DM interaction coefficient, D. Therefore, to 
estimate D in the continuum limit from the DFT calculation, 
we compute the long-wave length limit of the spin suscepti¬ 
bility, that is. 


Dp 


= 0 ) 

lim - ttt-r - 

g->0 lOqP 


( 2 ) 


Here, (a, / 3 ,7) = (x, y, z), {y, z, x), or (z, x, y) and Dp cor¬ 
responds to the coefficient for MadM-^/dfJ in Eq. ([T]i. Since 
we consider the skyrmions in the x-y plane under the total 
magnetic moment along the z-axis, hereafter we focus on Dx 
and Dy. In Eq. (O, we use the non-interacting spin suscepti- 
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bility defined as 


/,Z',Si ,S2,S3,S4 fe,m 

G'?sii'S2 *‘^™)^S2S3^i'S3iS4 {k + q,iujm+ii^i) ( 3 ) 

where, cr is the Pauli matrix and C?° is the non-interacting 
Green’s function in the orbital basis. Using this non¬ 
interacting spin susceptibility, we can write as Djs = 
i^/y)J2kDi3ik) with 


Dfiik) = lim --— 


<9 /(gn^fc+qr) ~ fi^nk) 

q—>-0 iOq^ gn^fc+q gnfc 

X {nfc|tT“|n'fc -I- q){n'k + q\a'^\nk), 


( 4 ) 


where |nfc) is the eigenvector of the Kohn-Sham Hamiltonian 
with the eigenvalue of e„fc. Hence, we can discuss the DM 
interaction in terms of the band structure. Although there 
is a sophisticated approach to compute IJ^, we employ the 
current simple approach to explore various parameters and 
materials. Furthermore, this approach is appropriate to obtain 
a guiding principle for controlling D as discussed below. 

Ab initio band structure. Figure [Ha) shows the DFT band 
structure of FeGe (black solid lines). Here, we include 
the spin-orbit couplings and assume the ferromagnetic 
moment along the z axis. The calculated local magnetic 
moment is 1.18 fiB per Fe atom, which is consistent with 
the experiments22i21 and previous calculations^. Using this 
electronic structure, we construct the tight-binding model 
made of Fe 3d and Ge 4p Wannier orbitals to reproduce the 
band structure below the Fermi level as shown in red broken 
lines. The densities of states for up spin (red line) and down 
spin (blue line) are also shown in Fig. [Tfa). As can be seen, 
there is a large exchange splitting, A. According to the energy 
difference of up and down spins for the Fe 3d orbitals, we 
obtain A = 1.17 eV. In Fig.[Tlb), the obtained tight-binding 
band structure around the Fermi level is illustrated with colors 
representing the weight of the up spin. Since we consider 
ferromagnetic electronic structure, each band can be basically 
characterized as either up-spin or down-spin band as shown 
in Fig.lHb). In addition, due to the spin-orbit couplings, there 
are several anticrossing points where complex spin texture 
emerges. 


Ab initio evaluation of the DM interaction. Let us start with 
the simp le two-band model in two dimensions considered in 
Ref.j^ The Hamiltonian is represented by a 2 x 2 matrix. 


H = kxcr^ + kya"^ + ma^. (5) 

We assume that the band dispersion is linear in the kx-ky 
plane, and dispersionless in the k^ direction. We introduce 
m to open a gap at the band crossing point as shown in Fig.|2] 
(a). The static spin susceptibility can be calculated analyti- 
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FIG. 1. Band structure of FeGe. (a) Comparison between DFT 
band structure (black solid lines) and tight-binding band structure 
(red broken lines). Densities of states for up spin (red line) and down 
spin (blue line) are also shown, (b) Detailed band structure around 
the Fermi level with colors representing the weight of up spin; that is, 
red (blue) lines correspond to up-spin (down-spin) bands. The Fermi 
level is set to zero. 


cally, and the result is 

Pi ^ — Ek) — S{fi + Ek) 

- n d{tqy) ~ AEk 

= -m) - e{-n - to)] (6) 

OTT 

where /i is the chemical potential and E^ = \Jkx + ky+ m?. 

It is interesting to note that Dy is negative for p, < —to, and 
positive for p > to as shown Fig. |2] (b). If we reverse the 
spin texture by modifying the Hamiltonian as ay —Cy, Dy 
becomes positive for p < —to, and negative for p > to. This 
result suggests that the position of the Fermi level and the spin 
texture around the anticrossing point are crucial to determine 
the sign of D. 

In real materials, the situation is not so simple as that of this 
two-band model. The anticrossing points form complex sur¬ 
faces in a four-dimensional space spanned by the energy and 
the wave number, and the electronic states around neighboring 
anticrossing points can hybridize with each other. In Fig. |2] 
(c), as a representative case, we schematically show how D 
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FIG. 2. Schematic pictures of the DM interaction around the 
hand anticrossing points, (a) Band structure of the two-band model 
defined in Eq. and (b) the chemical potential dependence of the 
DM interaction corresponding to this band structure. When two an¬ 
ticrossing points with different energies and spin textures with oppo¬ 
site chiralities reside close to each other, the peak structure appears 
as shown in (c). 



FIG. 3. Distribution of anticrossing points. Number of k points 
in 64x64x64 mesh where the up-spin weight, ui-|-, satisfies 0.4 < 
< 0 . 6 . 


changes as a function of the chemical potential when two an¬ 
ticrossing points with different energies and spin textures with 
opposite chiralities reside close to each other. Thus, when 
many anticrossings are densely clustered around the Fermi 
level, the sign and the size of D should change drastically. In 
fact, the band structure of FeGe has many anticrossing points 
around the Fermi level (see Fig. [T]i. To visualize the distri¬ 
bution of the anticrossing points as a function of energy, in 
Fig. 12 we plot the number of k points in 64 x 64 x 64 mesh 
where the up-spin weight, w^, satisfies 0.4 < w-f < 0.6. We 
can see that there are two peaks around /x = 0 eV and —0.5 
eV, where D is expected to change significantly. 

To analyze the contribution of each anticrossing point, next. 
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FIG. 4. Contribution of anticrossing points to the DM inter¬ 
action. (a) Band structure around the anticrossing points and k- 
dependence of Dx{k) at two different chemical potentials (b) p = 
—0.38 eV and (c) fj, — —0.34 eV at T = lOOK. Fermi surface is 
also shown in black lines. For comparison, the Berry curvatures at 
(d) fi = —0.38 eV and (e) p = —0.34 eV are plotted. 


we focus on the anticrossing points around the M (0 1/2 1/2) 
point. In Fig. 12 a), there are two anticrossing points between 
Y and M points, and between M and Z points. In fcj; = 0 
plane, such anticrossing points continuously form a closed 
loop around the M point. In Fig. |2b) and (c), we show the 
fc-dependence of Dx{k) when the chemical potential is be¬ 
low and above the anticrossing points shown in black lines in 
Fig. 12 a), respectively. We see that the texture around the M 
point drastically changes when the chemical potential sweeps 
across the anticrossing points, which is consistent with the 
simple two-band calculationa^i. Note that Dx{k) around the 
other Fermi surfaces such as the one between F-Y line are 
not negligible although the spin mixture is not so significant. 
This is because the effect of spin mixing extends away from 
anticrossing points; that is, D{k) has non-negligible values 
typically up to 0.2 - 0.3 eV away from anticrossing points. 
For comparison, in Fig.|4](d) and (e), we show the Berry cur¬ 
vature, which is the origin of intrinsic anomalous Hall 

conductivity (AHC)2ii2^. is defined as 




n^n' 


nh 


2hn{nk\vx\n' k) {n' k\vy\nk) 


( 7 ) 
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FIG. 5. Carrier density dependence of the DM interaction. DM 

interaction coefficients, Dx, and Dy, and the AHC, Oxy as a function 
of the carrier density, n. The inset shows the relation between the 
chemical potential /r and carrier density n. We use the rigid band 
approximation starting from the electronic structure of FeGe (n = 
0.0, /r = 0.0 eV). n = —1.0 (/r = —0.506 eV) corresponds to the 
carrier density of MnGe. 


where Vx and Vy are velocity operators. Since the spin mixing 
is important for both cases, there are common regions where 
the contributions to the DM interaction and AHC are large. 
However, in the Berry curvature, the summation in Eq. (|2]l is 
restricted Xo n' while it is not in Eq. (|4|i. As a result, only 
the restricted region is important for AHC, which is in sharp 
contrast to D. 

Eigure |5] shows the resulting Dx,y at T = 300K as a 
function of the carrier density together with the AHC, (Txy 
In the calculation, we use the rigid band approximation. The 
relation between the carrier density, n, and the chemical 
potential, p, is shown in the inset of Eig. |3 At /i = —0.506 
eV, the number of hole is 1.0 per Ee atom, which corresponds 
to the carrier density in MnGe. In Eig. |5l we find that Dx,y 
shows clear sign change from EeGe {Dx,y < 0) to MnGe 
[Dx^y > 0), which is consistent with the experimental sign 
change of skyrmion helicityi^, and is in sharp contrast to 
axy The positive (negative) hump structure in iT(/i) around 
/i ^ 0.4 (0.2) eV originates from the peak structure around 
/i ^ 0.5 (0.0) eV in Eig. |3 respectively. If we assume that 
contributions from p ^ 0(—0.5) eV is negative (positive), 
we can understand why D changes its sign around ji ~ 0.3 
eV. As for cjxy, the calculated value of axy in MnGe is larger 
than that in EeGe and there is no sign change. This behavior 
including size and sign agrees well with the experimental 
anomalous Hall contribution to axy for MnGe^ and for 

Mni_a:Eea:Ge22. 

Strain-induced huge anisotropy of the DM interaction. 

According to the above discussion, there are several ways 
to change the size and sign of D. As we have seen above, 
D can be efficiently controlled by carrier doping. We can 


also exploit the temperature dependence of the exchange 
splitting. When the relative position of up- and down- spin 
band changes, the distribution of anticrossing points in the 
band structure also changes, which will have a direct impact 
on D. As an example, we will show later the change in 
Dx,y of EeGe as a function of the moment per Ee atom in 
Fig- [11(b)- This mechanism can be related to the temperature 
dependence of the magnetic moment and skyrmion size in 
MnGe^^. Another interesting possibility is to make use of the 
strain effect. If we apply a strain to the system, the symmetry 
of the electronic structure can be lowered, and the distribution 
of the anticrossing points will change drastically. This effect 
is expected to be prominent especially when D changes its 
sign. Eigure |6] (a) shows the calculated Dx^y at n = —0.45 
for which we apply the uniaxial strain along the y direction. 
We find that the difference between Dx and Dy actually 
enhances particularly by the elongation along the y axis; Dx 
is about 40 % (400 %) larger than Dy at +2 % (h- 5%) strain. 
The carrier density dependence of the anisotropy, \Dy/Dx\, 
for fixed strain of +5% is shown in Eig.|6](b). We can see that 
the anisotropy becomes large particularly around the region 
of sign change. 

Discussion 

In the present calculation, we employed the rigid band ap¬ 
proximation. To examine its validity, we have performed a 
calculation for the other end material MnGe and doped neg¬ 
ative carriers by the rigid band approximation. As shown in 
Fig.lZl we have obtained a qualitatively similar result in that D 
is positive (negative) for the end material MnGe (EeGe). Thus 
the result that D for Mni_a;Eea;Ge changes its sign between 
a; = 0 and 1 should be robust, even when we go beyond the 
rigid band approximation. 

Regarding the crystal structure, it has been known that 
the magnetic moment for the optimized structure is much 
smaller than the experimental value within the local density 
approximation^^. On the other hand, for the experimental 
structure, the size of magnetic moment is similar to that in 
the experiment (~ l/ie)- In fact, within our calculations, the 
magnetic moment in EeGe decreases with decreasing the lat¬ 
tice constant as shown in Eig. 0(a), which is consistent with 
previous studies2122. As a result, D exhibits non-trivial mag¬ 
netic moment dependence as shown in Eig. 0(b). Since the 
magnetic moment dependence of D is also an important issue, 
let us next discuss how the change in the lattice structure or the 
magnetic moment affect D in Mni_a;Eea;Ge. In Eig. 0(c), we 
plot the lattice constant dependence of D, together with the 
distribution of band anti-crossings (d), (e) for lattice constants 
a = 0.99ao and I.OIoq, where oq is the experimental lattice 
constant. As can be seen, for larger a and magnetic moment, 
the energy difference of two peaks in the histogram of Eig. 0 
becomes larger (Eig.0e)). Consequently, the density at which 
D changes its sign (ric) becomes larger (Eig. 0c)). However, 
the qualitative feature of D is robust against the change in 
a. Therefore, our calculation successfully explains why D is 
positive (negative) for MnGe (EeGe), and provides a useful 
guideline for materials design of skyrmion crystal in metallic 
helimagnets. 
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FIG. 6. Strain dependence of the DM interaction, (a) Dx,y at 

Anisotropy of the DM interaction, \Dy/ \ at +5% strain along the y 
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I = —0.45 as a function of uniaxial strain along the y direction, (b) 
direction as a function of the carrier density. 
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FIG. 7. Carrier density dependence of the DM interaction start¬ 
ing from MnGe electronic strnctnre. Carrier density dependence 
of the DM interaction coefficients, and Dy. We use the elec¬ 
tronic structure of MnGe (n = 0.0) for the same crystal structure 
as that of FeGe and employ the rigid band approximation, n = 1.0 
(/i = 0.832 eV) corresponds to the carrier density of FeGe. 


TABLE I. List of parameters for FeGe and MnGe. Transition tem¬ 
peratures, Tjv (K), ferromagnetic exchange couplings, J (meVA^), 
helical periods, A(A), experimental DM interactions evaluated by 
D = Ttt J/A (meVA) and calculated D (meVA). J is evaluated by 
assuming J oc Tjv and the values of J = 52 meVA^ and Tjv = 30 
Kfor MnSi. 


Tn 

J 

A 

D (expt.) D (calc.) 

(K) (meVA^) (A) 

(meVA) 

(meVA) 

FeGe 278 

482 

700 

-8.7 

-10.1 

MnGe 170 

295 

30 

124 

107 


Finally, let us compare the quantitative values of D with 
experiments. For this purpose, we should rescale D using the 
exchange splitting as D = Using A = 1.17 eV for 

FeGe and A = 2.14 eV for MnGe based on the calculation 
shown in Fig.|7] we can obtain D = —10.1 meVA for FeGe 
and D — 107 meVA for MnGe. The experimental values 
of D can be estimated using J and the helical period, A as 
D = TttJ/A. Assuming that J scales to Tjv and using the 
values of J = 52 meVA^ and = 30 K for MnSU^, the 
experimental values of D for FeGe and MnGe are -8.7 meVA 
and 124 meVA, respectively, which are in good agreement 
with our results (see Table I). 

Method 

Crystal structure. In the calculations, experimental values 
are used for the crystallographic parameters^. For the 
pressure and strain calculations, internal coordinates of 
the atoms are fixed and only the lattice parameters are 
modified. To symmetrize and Dy, we use two different 
internal coordinates of atoms in the (4a) position, that is, 
{x,x,x) where ipe.Mn = 0.135 and a;Ge = 0.842, and 
the one with its 90-degree rotation along the z-axis and 
take the average. Note that our structure is right-handed 
according to Ref. 0. In the experiment, observed skyrmions 
on the right-handed crystal structure in FeGe (MnGe) are an¬ 
ticlockwise (clockwise), indicating that Dx,y < 0{Dy,,y > 0). 

DFT calculations. To evaluate Xo^ in FeGe, we perform 
the electronic structure calculation within the generalized- 
gradient approximation (GGA)^ based on the density 
functional theory^. We use ultrasoft pseudopotentials^ and 
plain-wave basis sets to describe the charge densities and 
wave functions with cutoff energies of 40Ry and 500Ry, 
respectively. We use 8x8x8 fc-point mesh. With including 
the spin-orbit couplings and assuming the ferromagnetic 
moment along the z axis, we obtain non-collinear magnetic 
structure. Using this electronic structure, we calculate Wan- 
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FIG. 8. Lattice constant dependence of the DM interaction, (a) Lattice constant dependence of the magnetic moment on each Fe atom for 
FeGe and (b) corresponding magnetic moment dependence of Dx,y (c) Dx,y with several lattice constants, a/ao = 0.98, 0.99,1.00, and 
1.01. Distribution of anticrossing points at (d) a/ao = 0.99, and (e) a/ao ~ 1.01. 


nier functions for Fe 3d and Ge 4p orbitals using wannierOO 
code^— . Based on the Wannier functions, we construct 
a tight-binding model on the restricted Hilbert space and 
calculate using 64 x 64 x 64 k-point mesh at T = 300K. 


The AHC is also calculated using the Wannier interpolation 
technique with 200 x 200 x 200 k-point mesh^. 
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